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Abstract. The recent determination of the dark matter density in the Universe by 
the WMAP satellite has brought new attention to the interplay of results from particle 
physics experiments at accelerators and from cosmology. In this paper we discuss 
the prospects for finding direct evidence for a candidate dark matter particle at the 
^ ' LHC and the measurements which would be crucial for testing its compatibility with 

C^- , cosmology data. 

o 

"^ ' 1. Introduction 

o 






Recent precision data on the Cosmic Microwave Background (CMB) and other 
astrophysical measurements have confirmed that a substantial fraction of the mass of the 



O-t. 

D . universe is in the form of dark matter; that is matter which which cannot be observed 

directly using conventional astrophysical techniques. The CORE |lj satellite provided 

evidence for intrinsic anisotropics in the CMB spectrum. The superior resolution of 

i^ ! the Wilkinson Microwave Anisotropy Probe (WMAP) mission has enabled the accurate 

5^ ! extraction of the baryonic and dark matter densities of the universe. This has confirmed 

that the dark matter density exceeds that of baryonic matter in the universe by a factor 

of six j2] . The mass density is expressed as a fraction (referred to as Q) of the amount 

required for an asymptotically fiat universe and it is presently known to an accuracy of 

7 %. The CMB results agree well with data obtained from supernovae and the rotation 

curves of spiral galaxies. 

Dark matter, as opposed to dark energy, can exist in two forms: non-luminous 
baryonic matter in the form of large planets or dead stars (MACHOS), or weakly 
interacting elementary particles that pervade large regions of space. Searches for 
MACHOS have concluded that these cannot be responsible for a substantial part of 
the dark matter 0. 

Particle physics candidates for dark matter fall into three basic categories depending 
upon their masses and interactions. Weakly interacting particles that are in thermal 
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equilibrium at a very early stage in the evolution of the universe will eventually fall out 
of equilibrium as the universe expands. The decoupling time (or temperature) when 
this occurs depends on the expansion rate of the universe as well as the couplings of 
these particles to other particles that are still in equilibrium. Particles that are (non- 
)relativistic at the time that galaxies start to form are referred to as (cold) hot dark 
matter. The simplest example of hot dark matter is a neutrino with a very small mass 
(< (9(20) eV) and of cold dark matter a very heavy neutrino (~ 100 GeV). In both cases 
the interaction rates are determined by the Standard Model of electroweak interactions 
(SM). 

The third type of particle dark matter can arise during the QCD phase transition 
as the universe cools down. In this case the result can be a gas of axions |1] . For certain 
values of the mass and coupling, these could account for a substantial part of the dark 
matter |2i . 

Hot and cold dark matter have implications for the distribution of visible matter 
(galaxies) throughout the universe jB] . The observed non-uniform distribution of such 
matter is believed to arise from the growth of quantum fluctuations remaining in the 
early universe after the end of inflation. These fluctuations are amplified by gravitational 
interactions and lead to the clumping that we observe today. The dark matter plays 
a vital role in this process as it dominates the total matter. Hot dark matter is free 
streaming and produces a universe with too little small scale structure. Cold dark 
matter is therefore strongly favored P|. 

There are models of particle physics that provide natural candidates for relic 
particles comprising the cold dark matter. This article will discuss these candidates 
and how they could be produced and their properties inferred from experiments at the 
LHC. We shall comment on the interplay between the LHC experiments and terrestrial 
experiments that aim to detect the dark matter directly. Most of our discussion focuses 
on supersymmetric candidates. 

Many of the candidate dark matter particles can be produced copiously at the LHC 
either directly or as decay products of the other particles. Measurements at the LHC will 
then provide information regarding the masses and couplings of both the dark matter 
candidate and the other particles with which it interacts and which are important in 
calculating its "relic" density. 

In section 12.11 we review methodology by which the relic density is computed and 
discuss how the dark matter density constrains the possible supersymmetric models. 
We then discuss the constrained MSSM (CMSSM) model for which most of the detailed 
simulations have been carried out. The consequences of relaxing some of the CMSSM 
constraints are then discussed. This section concludes with a discussion of dark matter 
candidates in gauge and anomaly mediated models and of the possibility that sneutrinos 
or gluinos could account for the dark matter. Section IHl then discusses the potentially 
relevant LHC measurements of supersymmetric particles. The possibility that particles 
arising from models of extra dimensions could account for dark matter is discussed in 
Section |3] Section |31 then summarizes the constraints that can be obtained from direct 
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searches for dark matter. In section ^ we illustrate how the measurements from LHC, 
astrophysics and direct detection complement each other and can be used to build up a 
complete picture of supersymmetric dark matter. 

2. Cold Dark Matter and Supersymmetry 

2.1. Supersymmetric Models 

Supersymmetry (SUSY) is a symmetry which relates fermions and bosons. In 
supersymmetric theories all existing particles are accompanied by partners having 
opposite spin-statistics. It is an appealing concept, for which there is currently no 
direct experimental evidence. As none of the partner particles have been observed, 
supersymmetry cannot be exact and must be broken so that the masses of the 
superpartners are larger than a few hundred GeV [Zj. Supersymmetry solves a long- 
standing problem with the SM, namely that the electroweak scale (~ the mass of the Z 
and W bosons) is of order 100 GeV even though radiative corrections in the model point 
to a scale that is many orders of magnitude larger, around the Planck scale (10^^ GeV ) 
where gravity becomes strong. A supersymmetric version of the SM is able to stabilize 
these radiative corrections and solve this "hierarchy problem" provided that the new 
particles have masses of less than a few TeV. 

When supersymmetric models of electroweak interactions are constructed, there 
often arises a symmetry, called R-parity, under which all the existing particles are even 
and the new super-partners are odd. If R-parity is violated then supersymmetric models 
can naturally give rise to unacceptably short nucleon hfetimes. If R-parity is conserved, 
there are two important consequences; the supersymmetric particles have to be produced 
in pairs and the lightest of them is absolutely stable. In supersymmetric models the 
Lightest Supersymmetric Particle (LSP) is invariably produced abundantly in the early 
universe and hence if R-parity is conserved these LSPs should still be present throughout 
the universe, providing a natural candidate for the dark matter. 

There are several candidates for the LSP, although many possibilities are ruled out 
by rather general considerations. Most importantly the LSP cannot carry electric charge 
as otherwise it would bind to ordinary matter, violating the very stringent limits which 
exist on exotic matter Pj. Remaining candidates are therefore the gravitino (partner 
of the graviton), G, the sneutrino (partner of the neutrino), z/, the gluino (partner of 
the gluon), g, and the lightest neutralino, x?; which is a mixture of the parters of the 
SM gauge bosons Z° and 7 and the Higgs boson. In this last case, the coupling of 
the LSP is determined by its content. It is most convenient to work in a basis where 
the content is regarded as Higgsino, wino and bino, the neutral fermion states of the 
supersymmetric Higgs sector and the SU{2) and U{1) supersymmetric gauge theories. 
The first component couples like a Higgs boson and therefore has very weak couplings 
to quarks and leptons of the first two generations. The second component couples more 
strongly than the third due to the greater SU{2) gauge coupling strength. 
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The simplest SUSY model, known as the "Minimal Supersymmetric Standard 
Model" or MSSM, is that containing one SUSY partner for each of the known SM 
particles together with two Higgs doublets to generate the masses of all the quarks 
and leptons. This model has over 100 parameters which control the masses and 
couplings of the new supersymmetric particles. These are determined by the mechanism 
responsible for SUSY breaking in the MSSM. There are several candidates for the 
breaking mechanism which produce qualitatively different mass spectra. In general it is 
assumed that there is an additional "hidden" sector of the theory with particles of very 
large mass where SUSY exists as a broken symmetry. This breaking is communicated 
to the partners of the SM particles via a mediation mechanism. 

Since gravity is known to exist, it is a natural candidate for the SUSY breaking 
mechanism and since it couples to all particles equally it is expected that all the sleptons, 
sneutrinos (partners of the charged leptons and neutrinos) and squarks (partners of 
quarks) will have a common mass at the energy scale where the mediation operates. 
Now the three couplings of the SM (the strong, weak and electromagnetic couplings) 
are precisely measured in current experiments, as are their evolutions with energy. If 
these couplings are extrapolated to very high energy assuming the existence of only the 
SM particles and their SUSY partners (with masses around 1 TeV) then it is found 
that they unify to a common value at a scale slightly below the Planck scale. This 
coincidence can be explained by postulating that the strong, weak and electromagnetic 
interactions are unified at this scale giving equal gluino and gaugino masses. These ideas 
are implemented in the SUGRA or constrained MSSM (CMSSM) model jHj. This has 
a total of five parameters: the common scalar mass niQ, which determines the squark, 
slepton, sneutrino and Higgs masses, the common gaugino mass ??ii/2, which determines 
the gaugino and gluino masses, tan/5 which controls the relative size of the couplings of 
the Higgs boson to up and down type quarks, the sign of the Higgsino mass parameter 
/i, and the common trilinear coupling Aq, whose value is relevant only in a few cases. 
It is a remarkable feature of this model that, as the parameters are evolved down to 
lower energies, the masses change in a calculable way and, due to the large coupling of 
the top quark to the Higgs, electroweak symmetry is broken spontaneously and masses 
for the W^ and Z'^ bosons generated automatically. As this model is simple, and fully 
described by a small number of parameters, it has been used most often for detailed 
phenomenological studies. Other models of SUSY breaking have different structures 
and we will comment on these later. 

2.2. Computing the Cold Dark Matter Density 

The maintenance of thermal equilibrium of particle densities in the early universe 
depends upon the interaction rates of the particles and the expansion rate of the universe. 
At very early times when the universe is at a very high temperature, all particles are 
in thermal equilibrium. Consider a particle that can annihilate against its anti-particle 
with a cross-section a. The products of this annihilation are other (lighter) particles. 
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The production and annihilation of the particles maintains their equilibrium. The time 
evolution of the number density n is described by: 

^ = -3Hn~ <av> {n' - n%) (1) 

The first term on the right represents the dilution caused by the expansion of the 
universe, with H being the Hubble constant. The second term represents the effect of 
the annihilation, with v being the velocity of the particles and < av > representing a 
thermal average of the cross-section times velocity. Ueq is the number density at thermal 
equilibrium given by a Boltzmann distribution. If a is sufficiently large, the second 
term ensures that n = n^q, and hence the relic density is exponentially suppressed at 
temperatures where the particles are non-relativistic. If a is very small then the particles 
cannot remain in equilibrium. For most processes interactions at sufficiently early time, 
or equivalently high temperature, are dominant, and equilibrium is maintained until the 
temperature falls below a freeze out temperature (Tp) when the interaction rate becomes 
too low and annihilation ceases. The number density at the present day is then given 
by neg{Tp) diluted by the Hubble expansion. The surviving number is most sensitive 
to a; a larger a ensures a lower number of survivors and a smaller contribution to the 
mass density Q. As the particle's mass increases its contribution to Q also increases and 
therefore a must be raised to satisfy experimental constraints. We shall frequently refer 
to this contribution to Q arising from a specific dark matter particle as Q^h'^, where the 
normalisation is appropriate for comparison with CMB data. 

Amazingly, if supersymmetric particles have masses of order 100 GeV as indicated 
by the fine tuning argument discussed above, their present abundance is of the right 
magnitude to account for the "observed" cold dark matter. It is this fact that has 
led to such a large focus on supersymmetric dark matter candidates. In practice, the 
computation of the relic density of LSPs is more complicated than indicated above. 
There can be several supersjTiimetric particles with comparable masses, all of which are 
falling out of thermal equilibrium at the same time. As well as particle-anti-particle 
annihilation, processes involving the scattering of different species are important ( "co- 
annihilation"). 

Several software packages have been developed to compute the cold dark matter 
density in a supersymmetric model with conserved R-parity ^01 E]- It is important to 
emphasize that the masses and couplings of the supersymmetric particles must be known 
before a prediction of the relic density can be obtained. Renormalization group equations 
(RGBs) are used to evolve the high scale CMSSM parameters and determine the soft 
parameters (masses and couplings) at the electroweak scale. Even small differences 
in these masses and couplings can result in significant changes in the predicted relic 
density. A recent survey found that the uncertainties in the masses of SUSY particles, 
defined by the differences between predictions from four independent programs, can be 
of order 20-30 % [T^. Roth DarkSusv [llj v. 4.0 and micrOMEGAs [10] v. 1.3 can 
be interfaced to several different codes for computing the full SUSY spectrum including 
ISASUSY [IH! and SUSPECT HH. 
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The main quantities which affect the rehc density are the mass of the LSP and its 
interactions, determined by its composition. Over much of the CMSSM parameter space 
the LSP is the hghtest neutrahno, which usually consists predominantly of the partner 
of the gauge boson of the ^(1) sector of the SM (the bino). The coupling constant 
of this is the smallest of the three SM couplings and the annihilation cross-section is 
therefore small. The dominant mechanism contributing to a is often XiXi ~~^ i~^i~ 
which includes slepton exchange and is sensitive to the mass of the lightest slepton. 
In the CMSSM therefore, as mo increases, a decreases and the relic density increases. 
Generically, supersymmetric particle masses cannot be too large if the universe is not 
to be overdosed by an excess of dark matter [Q > 1). 

In order to calculate the LSP relic density for all possible masses, one must 
consider a very large number of possible SUSY annihilation channels. In practice 
these contributions must be computed rapidly in order to explore wide regions of the 
multi-dimensional supersymmetric parameter space. An interesting solution, adopted 
by micrOMEGAs has been to rely on the CompHEP program [E! for computing the cross- 
sections. Loop corrections can be important in many cases, in particular those affecting 
widths of particles such as the Higgs. Use of one-loop results rather than a tree level 
treatment of the Higgs width can change the ^^^"^ ^^^^^ density value by a factor of 
two at large tan/3, due to the enhancement of the bb coupling. Thanks to efforts 
over more than a decade, the relic density can now be estimated from the electroweak 
scale MSSM parameters with an accuracy which in principle approaches 1%. However, 
supersymmetric threshold corrections will depend on the accuracy in determining the 
relevant SUSY particle masses. Also, the full QCD and the higher-order electroweak 
corrections need to be computed. If these are available by the time the LHC has collected 
a significant data set, the intrinsic accuracy on the dark matter density extracted from 
the supersymmetric particle spectrum should be comparable to that expected from the 
next generation of satellite-borne CMB experiments. 

As the SUSY particle masses increase, cosmologically viable solutions are possible 
only if the (co-) annihilation rates can be enhanced, for instance by accidental mass 
degeneracies. Furthermore, Xi annihilation through resonances may take place at large 
velocities, thus requiring a relativistic treatment. A formalism for relativistic thermal 
averaging was developed in Ref. [in] and has been adopted by several programs. 

If the LHC experiments are able to determine the masses and couplings of the SUSY 
particles, then the cold dark matter density can in principle be predicted using these 
software packages. A small change in the mass of the LSP may not have a significant 
effect on the LHC phenomenology but can have a dramatic effect on the cold dark 
matter density. Conversely, the masses of heavier SUSY particles such as squarks and 
gluinos are very important at the LHC, as these are usually the dominant source of 
produced SUSY particles, while their masses are usually irrelevant to calculation of the 
relic density. A consistent attempt to reconcile LHC data with that from dark matter 
observations will require that the same model can be applied to both. 
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Figure 1. Schematic plot of the mQ-mi/2 plane showing the regions where 
cosmologically interesting cold dark matter relic densities occur (light shaded regions) . 
The dark region at low ttiq is excluded by the requirement that the LSP be neutral 
(here it is a charged slepton), while the dark region at large mo is excluded by the 
requirement that electroweak symmetry breaking occur. If SUSY is to explain the 
possible anomaly in the measurement of the muon anomalous magnetic moment jl7| , 
the region to the left of the g — 2 line is preferred. A region at small mi/2 is excluded 
by the direct limits on the Higgs boson mass JHI- The contribution from SUSY to 
the observed ^Hl decay 6 ^ 57 excludes the region at small masses unless the SUSY 
contribution is cancelled by other new effects. 



2.3. Cold Dark Matter m the Constrained MSSM 

Within the constrained MSSM class of models the LSP is almost always the lightest 
neutralino. Cosmologically interesting LSP relic densities mostly occur only in certain 
regions of mQ-mi/2 parameter space illustrated schematically in Figure H The region 
at low mo and rrai/2 is referred to as the "Bulk Region". Here there are no accidental 
degeneracies of the SUSY particles, no enhancements of a and the LSP has a mass of 
less than 200 GeV. The region extending from the bulk region to large mi/2, running 
along the edge where the slepton becomes the LSP has an enhanced annihilation rate 
as the lightest slepton and LSP are almost mass degenerate. This is known as the "co- 
annihilation tail" . The region at large values of mo and mi/2 is sometimes known as 
the "rapid annihilation funnel" or "Higgs pole region" and occurs when the mass of the 
LSP is such that annihilation via an intermediate (s-channel) heavy Higgs boson (A) 
enhances a. Finally, at large mo there is the "focus point region" or "hyperbolic region" 
along the boundary beyond which electroweak symmetry breaking no longer occurs. 
We will now discuss in more detail the phenomenology of each of these regions. Note 
however that cosmologically interesting relic densities can also be generated in certain 
other scenarios, for instance when the LSP is nearly mass degenerate with the lighter 
stop squark at large Aq. 
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Figure 2. The mo-mi/2 plane for tan/3=:10 (left) and 50 (right) showing the region 
allowed by WMAP data and the existing constraints from accelerator data. The "L"- 
shaped dark (blue) and light (cyan) regions correspond to 0.094 < fl-^h"^ < 0.129 
(WMAP) and 0.10 < fi^/i^ < 0.30 respectively. The regions on the left of the vertical 
lines are excluded by the LEP-2 limits on the chargino (dashed) and the Higgs (dashed- 
dotted) masses. The region below the lower (brown) triangle is forbidden since here 
the LSP is charged. From Ref. [211 ■ 



2.3.1. Bulk Region and Co- Annihilation Tail The region at low values of mo and mi/2 
corresponds to the largest area of the CMSSM parameter space yielding a cosmologically 
interesting cold dark matter relic density. In this region, the SUSY particles are 
relatively light and therefore the region is severely constrained by negative results of 
searches for super symmetric particles at LEP-2 and the Tevatron. In the CMSSM, the 
Higgs mass is constrained to be less than 134 GeV [211! • Searches for the Higgs boson at 
LEP result in a lower limit on the Higgs mass of 114.4 GeV JH]- This and the lower limit 
on the chargino mass [Zj reduce significantly the allowed portion of this bulk region. 
Virtual effects from SUSY particles can affect the rates for certain rare processes. In 
particular, the observation of the loop-mediated 6 — >■ 57 process provides a constraint 
which excludes small masses and is particularly important at negative /i and large values 
of tan /3. If SUSY is to explain the possible anomaly in the measurement of the muon 
anomalous magnetic moment ^\, small masses are preferred. 

As the masses of the SUSY particles increase, acceptable relic density is obtained 
by enhancing the annihilation rate. This happens in a narrow strip in the mo-mi/2 
plane (the co-annihilation tail), where the x^f -^ t^ co-annihilation process becomes 
enhanced due to near degeneracy between the masses of the f and x? LSP. The precision 
of the WMAP measurement of i^^^^ ^^^ reduced the cosmologically allowed region 
to a very narrow one for each value of tan/3 and sign of /x. The upper tip of this 
region, corresponding to the heaviest supersymmetric particle spectrum, is defined by 
its intercept with the slepton LSP (m^ < m^) boundary. This is typically located at 
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Figure 3. The focus point region in the CMSSM ttiq — rni/2 plane for Aq = 
0,tan/3 = 10, fj, > 0. The medium (green) region is consistent with WMAP rehc 
density constraints (0.094 < fi^/i^ < 0.129), while the light (yellow) region is excluded 
by LEP2 bounds. From Ref. [21. See also Ref. J2S1- 



mi/2 ~ 0{1 TeV). As we shall see below, this region can be entirely covered by LHC 
searches. 

2.3.2. Focus Point Region In the focus point region at high m^ the LSP changes its 
nature. Instead of being dominantly gaugino, as is the case elsewhere, it acquires a 
significant Higgsino content as /i is driven to small values j22] • At high values of tan j3 
the region where this occurs can lie at rather low values of mi/2- As the Higgsino 
component can have a coupling to the SM gauge bosons, which is forbidden for the bino 
component by gauge invariance, a is enhanced |23j . Furthermore the x? becomes nearly 
mass degenerate with the Xi and X2- Consequently a number of additional annihilation 
and co-annihilation processes contribute which together can reduce the x? relic density 
further to levels compatible with CMB data. Unlike the co-annihilation tail, the focus 
point region can extend to very large values of m^ leading to squarks, sleptons and 
gluinos that are too heavy to be observed at LHC. 

The very existence of a cosmologically acceptable focus point region in CMSSM 
parameter space is sensitively dependent on the parameters describing the model, 
especially the top quark mass. In addition the large value of rriQ results in a "fine- 
tuning" problem; a delicate cancellation is required between different contributions to 
the Higgs potential to ensure that Mz « niQ. This problem worsens as mo increases. 



2. 3. 3. Higgs Pole Region When the value of tan j3 exceeds ~ 30 another opportunity for 
enhanced Xi annihilation is offered by the presence of heavy neutral Higgs particles, A'^ 
and H'^ with masses m^ such that m^o ~ ?7i^/2. In this case the annihilation of lightest 
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Figure 4. Xi LSP relic density in the mo — mi/2 plane for tan/3 — 5 and differing 
gaugino mass assumptions. The cosmologically allowed relic density regions (light 
(yellow and green) shading) are displayed for r — I (left panel) and for r — 0.6 (right 
panel). Dark (brown) shaded regions in the lower right are excluded due to a charged 
LSP. Medium (pink) shaded regions on the left are excluded by LEP bounds on the 
chargino mass. The Higgs mass contour of nih = 113 GeV is given by the dotted 
vertical line. Adapted from Ref. 27 ) 



neutralinos is enhanced through resonant (s-channel) heavy Higgs exchange [^ . These 
processes can allow CMSSM models with mo and/or mi/2 of order 1 TeV or greater to 
satisfy relic density constraints. The cosmologically acceptable region in these models 
takes the form of a "funnel" pointing toward large tua defined by relic density contours 
passing around the annihilation pole. For values of tan/?^60 the coupling of the Higgs 
boson to bottom quarks becomes large and reliable computations are not possible. 



2.4. Beyond the CMSSM 

As the assumptions that lie behind the SUSY breaking in the CMSSM are relaxed, 
the range of SUSY particle masses for which there exists an acceptable dark matter 
candidate is considerably expanded. Within the CMSSM scheme, the assumption of the 
unification of the gluino and gaugino masses at a high scale ensures that the neutralino 
LSP is predominantly bino in the bulk region. If the LSP were dominantly wino on 
the other hand, it would have larger couplings and larger annihilation rates. Its mass 
could therefore be larger. This is illustrated in Figured from Ref. [2Zj which shows the 
effect of changing r = M2/M1 from its unification value of r = 1. For r = 0.6 the LSP 
becomes wino-like and the bulk region expands as shown in the right hand plot. 

The CMSSM also assumes that the soft scalar masses have a common value at 
the unification scale. This assumption can also be relaxed leading to modified relations 
between the squark, slepton and Higgs boson masses. This causes the LEP Higgs mass 
limit to no longer directly constrain the slepton and squark masses. The bulk region is 
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not greatly affected by this modification, and the co-annihilation tail still exists, however 
the position and character of the focus point region can be radically altered since the 
squark and slepton masses no longer have to be large to obtain a small value of /i |28j . 

Another mechanism of SUSY breaking, "anomaly mediation" (AMSB) j31] , results 
in a SUSY spectrum where the LSP is a wino. The simplest version of this class of models 
is ruled out as it predicts tachyonic sleptons, however variants can escape this constraint 
P3] . The characteristic feature of AMSB models is the near degeneracy of the LSP and 
the charged wino resulting in a large annihilation cross- sect ion. Consequently if the LSP 
has a mass around 100 GeV far too few survive today to make a significant contribution 
to the dark matter. Nevertheless if the mass of the LSP is very large it is possible that 
a satisfactory dark matter density could be achieved. 

In addition to gravity mediated models such as the CMSSM there also exists 
a class of gauge mediated SUSY breaking (GMSB) models jSHl IHOj in which the 
supersymmetry breaking is mediated by gauge interactions. GMSB models assume that 
supersymmetry is broken with a scale a/F in a sector of the theory which contains heavy 
non-SM particles. This sector then couples to a set of particles with SM interactions, 
called messengers, which have a mass of order M. The mass splitting between the 
superpartners in the messenger multiplets is controlled by a/F. One (two) loop graphs 
involving these messenger fields then give mass to the superpartners of the gauge bosons 
(quarks and leptons) of the SM. GMSB models are preferred by some because the 
superpartners of the Standard Model particles get their masses via gauge interactions, so 
there are no flavor changing neutral currents, which can be problematic in the CMSSM. 
The non-observation of supersymmetric particles and new particles associated with the 
messenger sector provides a bound on the SUSY breaking scale: F^^m-+M GeV^ 
where 02 ~ 1/28 is the coupling strength of the SU{2) electroweak interactions. Using 
M, m-+ > 100 GeV implies vF^2 TeV however in most models F is at least ten times 
larger than this. Note that for gauge mediation to dominate F cannot be too large. 

In GMSB models the characteristic spectra of superparticles are different from those 
found in CMSSM models. In particular the lightest supersymmetric particle is now the 
gravitino (G) whose mass is given by uIq ~ — ^\r ^V, implying that ttiq > 0.2eV. The 
gravitino can have a mass of up to a GeV or so in specific models. It has feeble couplings 
and can be produced with significant rates only in the decays of particles which have 
no other decay channels. Its interactions are so weak that annihilation in the early 
universe typically only occurs via gravitational interactions and it falls very rapidly out 
of thermal equilibrium |31| leading to overclosure of the universe. This problem can be 
avoided if the temperature to which the universe is reheated after the end of inflation is 
sufficiently low to avoid over-production P2]. An alternative source of gravitino dark 
matter is non-thermal production in which the next to lightest SUSY particle (NLSP), 
which is usually a slepton or bino, decays to the gravitino. The lifetime of the NLSP can 
be sufficiently long for it to be effectively stable as the universe evolves, decaying only 
later to produce gravitino dark matter, long after it has fallen out of thermal equilibrium 
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In certain SUSY models it is possible for the LSP to be a sneutrino. This is however 
not favored in most models as the sneutrinos are usually heavier than the right handed 
sleptons due to the their larger gauge couphngs. Nevertheless a stable sneutrino is a 
good dark matter candidate fH^ providing its mass is around 500 GeV ^7\. Interaction 
rates of sneutrinos in direct detection experiments (see Section Ej) are expected to be 
similar to those of heavy neutrinos of the same mass. Consequently current bounds 
from these experiments rule out sneutrino dark matter for sneutrino masses less than 1 
TeV ISH]. 

There are a small number of SUSY models in which the gluino is the LSP [39j. 
however these are rare since ordinarily the larger coupling strength of QCD drives the 
gluino to larger masses than the other gauginos. Gluino LSPs produced thermally 
in the early universe are also unlikely to contribute significantly to the dark matter 
because they annihilate very efficiently via strong interactions leading to a number 
density typically of order 10~^° that of baryons. Gluinos cannot therefore account for a 
significant part of the dark matter unless they are produced non-thermally. 

So far we have mostly assumed that the LSP was initially in thermal equilibrium. 
If this was not the case, then the prediction of its relic abundance is no longer valid. In 
particular it has been pointed out that the decay of "moduli fields" expected in some 
unification models can produce too many LSPs pi)]. In gauge mediated scenarios, the 
excess production of gravitinos restricts the range of allowed masses pT|. In the case of 
anomaly mediation, it has been suggested that this mechanism can produce sufficient 
LSPs to obtain the correct relic density for LSP masses around 100 GeV 02] • 

3. LHC measurements of SUSY Dark Matter Properties 

3.1. Inclusive Searches 

The first issue is the reach of the LHC experiments for discovering supersymmetry. The 
reach is dependent on the details of the SUSY model, however the characteristic signals 
of missing transverse energy (from invisible LSPs), a large multiplicity of hadronic jets 
(from decays of squarks and gluinos) and/or isolated leptons (from decays of sleptons 
and gauginos) are sufficient to ensure detection over a large mass range. The left plot 
in Figure El shows the expected sensitivity of the CMS detector via the Jets + EJp'^^'^ 
search channel to models in the CMSSM mQ-mi/2 plane as a function of the integrated 
luminosity. About 30 fb~^ of data ensure that a signal can be observed over the whole 
extension of the co-annihilation tail. The right plot of Figure El shows the sensitivity 
with 100 fb~^, in a variety of channels characterised by the number of observed isolated 
leptons. The most sensitive channel is that with jets and missing transverse energy and 
no isolated leptons labelled "01" . Signals with two isolated leptons of the same charge 
(marked "21 SS") can arise from the pair production of gluinos and their subsequent 
decay to charginos. Since the gluino is a Majorana fermion, it can decay into charginos 
of either charge. Their decays then give rise to this signal. Signals with two isolated 
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Figure 5. Sensitivity of LHC experiments to supersymmetry in the CMSSM mo-mi/2 
plane. Left: reach of the Jets + i?™'"^ search channel to obtain a 5 cr sensitivity 
for different integrated luminosities shown as solid lines. The dashed lines show fixed 
values of the squark and gluino masses. Right: reach of different search channels for 
100 fb^^ of data indicated by the solid lines The cosmologically preferred region is 
indicated by the dark (blue) line extending over the excluded region where the LSP is 
charged (adapted from Ref. 03]). 



leptons of opposite charge (marked "21 OS") arise from the same source but also from 
decays of the type X2 ~^ ^^Xv Trilepton events (marked "31") arise predominantly from 
events containing X2 ^^^d ^ chargino. 



3.2. Measurement of Cold Dark Matter Properties 

The observation of signals characterized by large missing transverse energy, jets and/or 
isolated leptons at the LHC will give a strong indication that supersymmetry or some 
other New Physics resembling SUSY such as UED (Section 0)) is present. It will also 
provide circumstantial evidence that R parity is conserved. The limited sensitivity of 
LHC experiments to the LSP lifetime {t^ > 1 ms ^3]) means that final confirmation 
that cold dark matter is supersymmetric in nature will require observation of signals 
in astroparticle experiments consistent with mass and cross-section predictions derived 
from LHC measurements. This is especially important given that the dark matter 
relic density measured by CMB experiments could consist of non-supersymmetric 
components in addition to that provided by the LSP. It will therefore be vital to measure 
sufficient SUSY parameters to verify consistency of the observed signals not only with 
CMB data but also with signals in these non-accelerator experiments. 

Initially these measurements will be performed within the context of a specific 
SUSY breaking model, such as the CMSSM, by comparing significances in the various 
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inclusive Jets + EJp^'^^ + n leptons channels discussed above. As more data accumulates, 
it will be possible to begin to reconstruct SUSY particle decays. This may enable 
model-independent measurement of the masses of some SUSY particles. However this 
will not necessarily be sufficient to provide a model-independent prediction of the dark 
matter relic density or other dark matter properties. This is because we do not know a 
priori which annihilation and co-annihilation processes dominated in the early universe, 
or equivalently in which region of which parameter space Nature's chosen model lies. 
With only limited knowledge of the SUSY mass spectrum it will therefore be necessary 
initially to assume a specific SUSY breaking model. This will enable the full SUSY 
mass spectrum to be obtained from a limited set of measurements, with additional 
input coming from the unification assumptions. Only once detailed model-independent 
measurements of more of the individual SUSY masses and couplings have been made 
can this model-dependency in the estimates of dark matter properties be reduced or 
eliminated. 

As SUSY events will each contain two invisible LSP's, it will not possible to fully 
reconstruct the final states in these events. Measurements using exclusive channels will 
therefore focus on identifying key features of the SUSY particle decays. This typically 
involves measuring the end-points and shapes of invariant mass distributions of leptons 
and jets arising in the decay of SUSY states j^ ^1 EHj- Measurement of CMSSM 
model parameters can be accomplished through a global fit of the parameters to the 
positions of these edges, while model-independent mass measurements can be obtained 
by solving the various mass relations simultaneously. 

The exclusive measurements which can be used to estimate dark matter properties 
in either a mo del- dependent or model-independent fashion depend strongly on the (co- 
) annihilation processes which drove the reduction of the dark matter density in the 
early universe. Within the context of the CMSSM the regions of parameter space where 
these different processes dominate correspond those discussed in Section ITHl We shall 
therefore now discuss the measurements which can be made in each of these regions. 

3.2.1. Bulk Region and Co- Annihilation Tail In the bulk and co- annihilation regions 
of CMSSM parameter space the mass of at least one of the sleptons or staus lies 
between that of the X2 ^^'^ ^^^ Xi (LSP). The production of supersymmetric particles 
is dominated by squarks and gluinos as these have strong interaction couplings. If 
gluinos are heavier than squarks then they decay to squarks. In these regions there is 
a significant decay rate for q — > qx2 ~^ qi'^i^ -^ Q'^^^^Xi and/or the equivalent stau 
channel. This decay chain gives rise to a jet, from the quark, a dilepton pair of the same 
fiavor and opposite charge, and missing energy, from the Xi- After events are selected to 
contain a pair of isolated leptons, jets and missing energy, the left plot in Figure El shows 
the dilepton invariant mass distribution for a typical model. A clean structure is visible 
with a kinematic end-point arising from this decay chain. Events beyond the end-point 
are due to other SUSY particle decays. The position of this end-point is determined 
by the masses of £^, X2 ^^^ Xi- If events are selected with a dilepton pair below this 
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Figure 6. Invariant mass distributions obtained from the LHCC Point 5 benchmark 
model (too = 100 GeV, TOi/2 = 300 GeV, Aq = 300 GeV, tan/3 = 2.1, fi > 0) for 
100 fb""'^ of data. Left: the invariant mass of a pair of opposite sign same flavor 
leptons arising from squark decay. Right: the invariant mass of a pair of opposite sign 
same flavor leptons and jet arising from squark decay. From Ref. j44| . 



end-point and this pair is then combined with all possible jets in the event, the mass 
distribution of the jet-dilepton combination with the lowest mass displays a kinematic 
structure whose end-point depends also on the quark mass (Figure IHl (right)). There is 
also a minimum value of the mass of the dilepton and jet and an upper value of the 
mass of a jet and single lepton. 

These four quantities can be used to solve for the four masses: q, i^, X2 ^^^ Xi- 
It is important to note that a model is not required; simply the decay chain must 
be identified. When the masses are extracted precisions ranging from ~ 3 % for the 
ql to < 14 % for the Xi fo^^ the LHCC Point 5 benchmark model are obtained |47j . 
Note however that the errors are strongly correlated as is illustrated in Figure [7| an 
independent measurement of one mass improves the precision on all of the others. Note 
also that when stau co-annihilation dominates, the above signatures can involve very soft 
taus which can be difficult to observe or measure with significant accuracy. Nevertheless 
if the mass difference between the fi and either of the lighter neutralinos is more than 
a few GeV then it should be possible to make at least a rough estimate of the relevant 
end-point positions. 

Initial model-dependent (see Section ing predictions for dark matter properties can 
be obtained from measurements such as these in two different ways. One approach is 
to fit the parameters of the assumed SUSY breaking model to the measured model- 
independent SUSY particle masses. The model parameters thus obtained can then be 
varied within their statistical and systematic errors to produce a probabihty density 
function for e.g. Q^h^ (left hand plot in Figure IH}- This approach however does not 
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Figure 7. Plot showing errors on masses determined from tlie distributions in FigureEl 
Note the strong correlation. From Ref. |44j . 



take into account the correlations between the errors on the different measurements. 
An alternative approach which does take these correlations into account is to perform a 
global fit of the model parameters directly to all the end-point measurements prior 
to calculating the Vt^h?' PDF (right hand plot in Figure IH}- For the SPSla bulk 
region benchmark model J3H] a recent study [1^ has shown that with 300 fb^^ of 
data statistical measurement precisions ~ 2%, 0.6%, 9% and 16% can be obtained 
for the CMSSM parameters mo, ?tii/2; tan/5 and Aq respectively. The sign of /i is 
also well-determined from the fit. It is estimated that Vt^h"^ can be pedicted with a 
statistical(systematic) precision ~ 2.8% (3.0%). Further dark matter properties such as 
the spin-independent Xi-nucleon elastic scattering cross-section agi (see Section EJ are 
estimated with statistical precisions ^1%, far smaller than the estimated factor > 2 
systematic uncertainties in these quantities. 

Full model-independent estimation of the dark matter density will require use of 
the model-independent mass values discussed above together with other measurements 
of quantities such as the masses of the lighter staus, heavy Higgs bosons and heavier 
neutralinos in order to show that slepton or stau (co-) annihilation dominates in the 
early universe (i.e. that the model does not lie near the heavy Higgs pole, or in the 
Focus Point region). More work on this subject is required. 



3.2.2. Focus Point Region The phenomenology of the focus point region differs from 
that of the bulk and co-annihilation regions primarily due to the large masses of 
the sfermions. Depending on the parameters, the LHC may perform measurements 
providing some constraints on the dark matter density, or only observe an excess of 
events due to supersymmetry, or even, at the upper end of the parameter range, not 
observe supersymmetry at all. Sequential two-body, slepton-mediated, decays of X2 to 
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Figure 8. Probability density functions for the dark matter density as can be 
determined at the LHC. The left-hand distribution shows dark matter densities 
calculated from post-LEP benchmark model B' [HOI by first fitting to invariant 
mass edges with model-independent SUSY particle masses, and then fitting CMSSM 
parameters to the measured mass distributions, before using those parameter values 
to calculate D.^h'^. The range 0.094 < fl^h"^ < 0.129, favoured by WMAP, is shown 
for comparison. The right-hand distribution was obtained from the SPSla benchmark 
model j42j by directly fitting CMSSM parameters to the invariant mass edges before 
using those parameter values to calculate il^h^ (from Ref. (49) ) . 



Xi, which form the basis of the measurements discussed previously, are kinematically 
suppressed. Some other signatures may still observed in these models, depending on the 
value of Am = m^o — m^o . 

A2 Al 

In models with Am < rriz, such as the LHCC Point 4 model (mo = 800 GeV, mi/2 
= 200 GeV, tan/? = 10, /i > 0, Aq = 0; not strictly speaking a focus point model but 
a good example of a model with niQ » mi/2) studied in Ref. [SHIll], the X2 decays 
almost exclusively to three-body final states involving a Xi and two opposite-sign same- 
flavor leptons or quarks. The dilepton signature is easiest to identify and provides 
an invariant mass end-point sensitive to Am which can be measured with flavour 
subtraction techniques 01]. Combination of the dilepton pair with a hard jet may 
also give sensitivity to the ql mass while rriQ is sufficiently small for q^ production to be 
significant. The mixed Higgsino-gaugino nature of the heavier neutralinos and charginos 
leads to significant couplings to the Z^, and h^ bosons, giving a characteristic Z^ peak in 
the dilepton mass spectrum. The production cross-section and px distribution of these 
Z^ bosons can be used to constrain heavier neutralino and chargino masses jUj. Further 
mass constraints were also obtained from a dijet invariant mass end-point sensitive to 
rrig — rn^o, albeit with potentially large combinatorics. Sensitivity to squark and other 
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sfermion masses is however difficult to obtain. 

Model-dependent techniques, such as the analysis of the shape of the E^"^ 
distribution of SUSY events, are then required. A combined fit of CMSSM parameters 
to the results of all the available measurements for LHCC Point 4, gave a precision 
~ 6% for mo, ~ 1% for mi/2 and ~ 15% for tan/? (300 fb~^). Unsurprisingly, the 
rriQ measurement precision was relatively poor compared with that for m,i/2 due to the 
difficulty of identifying and measuring the sfermion signals. The lack of sensitivity to 
Aq and sign(/i) in this study indicates that only an order of magnitude estimate of Q^h"^ 
could be obtained. 

In the parameter region where Am > mz two-body decays of X2 ^o Xi + Z'^ ^-^id, 
for Am > nifi, Xi + h become accessible. Phase-space suppression of the three-body 
X2 modes and dominant X2 production in complex g three-body decays involving heavy 
flavors means that, even if an inclusive signal is observed, it is difficult to find mass 
distributions providing measurable end-points with little combinatorial background. 
Parameter measurements must therefore rely upon other techniques involving the Z° 
or h signatures. This is an area where further study is still required to assess the LHC 
potential. 

The large masses of the strongly interacting particle partners cause an increasing 
fraction of the total SUSY production cross-section to be taken up by direct production 
of neutralinos and charginos. These events are characterised by little jet activity with 
instead multi-leptons, or soft b-jets, arising for instance from Z^ or h decay. Interest 
in specific searches for these signatures has been re-awakened by results of a recent 
study j221- It indicates that the reach in the focus point region of gaugino searches at 
a 1 TeV e+e" linear collider, where gaugino pair production can be detected easily, can 
be greater than that of the LHC, when using the conventional jets -|- EJp^'"^^ channels. At 
the LHC the production cross-sections for the electroweak gauginos in the focus point 
region are comparatively large, and the SM backgrounds, with appropriate cuts, can be 
made much smaller. The disadvantage of these channels is that hard cuts on jet activity 
are required to reject hadronic backgrounds leading to low efficiency and potentially 
large systematics due to poor understanding of mini-jet backgrounds from pile-up and 
the underlying event. These searches resemble those performed at the Tevatron. The 
characteristic signatures are dileptons or trileptons -|- E'™'^* arising from the production 
of electroweak gauge bosons in gaugino decay. In an early study of these channels ISHl 
a plausible central jet pr cut of 25 GeV was assumed, consistent with that used in more 
recent central jet veto studies for Higgs searches in the vector boson fusion channel |54j . 
Only limited sensitivity to the high mo region was found. However with a more refined 
analysis and a finer granularity scan of the narrow focus point region it may be possible 
to improve the reach of the method. 

At very large values of mo and mi/2, the Xi becomes almost pure Higgsino and 
the xf can be almost mass degenerate with it (the so-called "inversion region" of the 
hyperbolic branch .25]). In this scenario, m-± — m^o and m^o — m^o are once again 
small and the two-body gauge and Higgs boson decay channels are closed. Three-body 
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decays again dominate, however the chargino and neutrahno mass differences can be 
such that only very low pt leptons and jets are produced (although the NLSP does 
decay inside the detector |2S1)- Triggering upon such events is likely to be very difficult 
at the LHC, and measurement of sparticle masses still more so. The limited scope for 
SUSY particle mass measurement in this region means that estimation of the Xi I'^lic 
density is a difficult task. This problem is compounded by the fact that in this scenario 
the value of Q^h"^ is strongly dependent on the top mass, which will only be known to 
a systematics-limited accuracy ~ 2 GeV at the LHC, after one year of low luminosity 
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Following any SUSY discovery at the LHC, the first indications that the focus point 
scenario is realised would likely come from a comparison of signal significance in different 
inclusive channels distinguished by the required number of leptons. In particular the 
signal significances in the jets + EJp'''^^ + n > 1 lepton channels are enhanced relative 
to that in the n = leptons channel in the focus point region due to the abundant 
production of heavy flavors in gluino decay and leptons in chargino and neutralino decays 
via gauge bosons. Further evidence could then be provided by exclusive channels. At 
low Am evidence for a di-lepton/di-tau end-point can be observed with a lepton/tau pt 
asymmetry consistent with three-body decays jlS|; while at high Am no evidence for 
such an end-point would be observable. 

Relic density estimation within the CMSSM framework at low Am would require 
use of the mass measurements and model parameter fitting techniques discussed above. 
At high Am, mi/2, tan/3 and sign(/x) could potentially be estimated from measurements 
involving direct chargino/neutralino production. But without an accurate estimate of 
mo, the magnitude of f2^/i^ would be difficult to assess. 

Estimation of the Xi relic density in a more model-independent manner in the 
focus point scenario is even more difficult. Quantities which would need to be measured 
include the neutralino mass matrix parameters Mi, M2, tan/3 and /i, and tua- While 
it may be possible to measure some or all of the neutralino/chargino sector parameters 
from direct production studies (high Am) , the A may well be inaccessible to direct study 
(via e.g. H/A -^ tt) due to its very high mass at high m-Q. However, if H/A — > tt 
were observed then this would enable a rather accurate measurement of tan/5 to be 
performed [HE] • 

3.2.3. Higgs Pole Region The potentially very large masses of the strongly interacting 
super-partners in these "rapid-annihilation" models means that, if realized in nature, 
their discovery at the LHC is by no means guaranteed (see, for example, the benchmark 
model "M" in Ref. [SOD- Nevertheless, if their spectrum is within reach of the LHC, 
it may be possible to measure many of the model parameters through the use of 
invariant mass end-point techniques. The rapid-annihilation funnel region in general 
occurs in the high tan/? mo — mi/2 plane relatively close (in all but the highest tan/? 
scenarios) to the stau co-annihilation tail discussed in Section K12.il Consequently for 
many rapid-annihilation models the fi is light and nif-^ < m^o. Ditau signatures may 
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therefore be used to identify SUSY events containing Xg decays, and fits to invariant 
mass thresholds and end-points involving one or more of these taus used to constrain 
either CMSSM parameters or model-independent SUSY particle masses. This analysis 
should be somewhat easier than in the stau co-annihilation case (for the same total 
SUSY cross-section) since here the fi — x^ mass difference is greater and hence the taus 
from the fi decays have greater pt- 

In addition to ditau production further signatures can be used to identify and 
measure SUSY events in rapid-annihilation models. These models in general possess 
relatively large values of ttiq compared with those of models in the stau co-annihilation 
tail, and this enhances the Higgsino content of the lighter neutralinos and hence the 
branching ratio of X2 ~^ Z^/h + Xv This allows further constraints to be placed on 
masses through identification of Z^ and h^ production in high E^^^^ events jSZj and then 
measurement of invariant mass end-points in Z^q and hq invariant mass distributions. 
Additional constraints can be obtained by searching for g^ -^ XiQ decays producing 
dijets -|- E^^^^, which give sensitivity to the qn mass |23EZI- Given an estimate of the 
mass of the Xi from e.g. the ditau analysis the measurement can be performed either 
by analyzing the shape of the pr distribution of each of the jets [SIl Elj or by fitting 
to the end-point of the distribution of the dijet mT2 variable jUl EHl- 

A case study of CMSSM parameter fitting in a model with this phenomenology is 
the analysis of LHCC Point 6 (mo = 200 GeV, my2 = 200 GeV, Aq = 0, tan/3 = 45, 
/i < 0) described in Refs. [Snil^. By using a combination of measurements of ruh and 
rriq^, constraints on rrig and rrig — mi obtained from g decays to 616, and fits to the ditau 
end-point, measurement precisions ~ 12%, 5% and 4% were obtained for mo, ''^1/2 and 
tan/3 respectively. However, the analysis was relatively insensitive to Aq and the sign 
of /x, which would limit the predictivity of these data in estimating the relic density. 

In models with very large tan/3 it is possible for the rapid- annihilation funnel to be 
sufficiently wide for the cosmologically-allowed region on the high m^ side of the base 
of the funnel to lie in the region where nif^ > m^o. The phenomenology of such models 
resembles more closely that of the focus point models discussed in Section IH.2.2| with 
characteristic Z^/h + E^'^^'^ {Am > rriz or rrih) or ditau + E^'^^''^ (Am < mz) signatures 
arising from X2 decay. 

Given a set of measurements of CMSSM parameters from the above analyses it 
should be relatively straightforward to calculate the Xi relic density, although the 
resulting error could well be significant due to its strong dependence on both SUSY 
and SM {rritop) input parameters |2i|. In order to obtain a more model- independent 
estimate of the relic density it will be necessary first to confirm that Xi annihilation in 
the early universe was dominated by resonant s-channel A production. Assuming that 
the x^ mass has been determined through the ditau -|- i?^**'^ analysis m^ must now be 
measured independently through a search for e.g. H/A -^ tt or (better) H/A — *• fifi 
to check that rriA — 2m^o. It is interesting to observe that in the CMSSM parameter 
region where a cosmologically interesting CDM density is achieved due to the rapid 
XX annihilation process through the A^ pole, the LHC should observe the A'^ and H^ 
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Figure 9. Estimated number of Higgs bosons observable at the LHC in the Ma, 
tan/3 plane. The cosmologically interesting solutions belonging to the A funnel regions 
obtained with a parameter scan are shown by the dots. In the relevant region of 
parameters, the LHC will be essentially observing all the Higgs bosons (modified 
from (61j ) . 



bosons for masses at least as high as 800 GeV (60], corresponding to mi/2 ~ 950 GeV 
with 300 fb~^ integrated luminosity (see Figure IH)). An LHC upgrade, increasing the 
luminosity by an order of magnitude, would ensure that the whole region is covered. 
The A mass measurement accuracy depends crucially on the availability of the A° -^ /z/z 
channel in which case a (9(0.1%) precision can be reached with 300 fb~^ {5m/ m ~0.1% 
for tua = 300 GeV at tan/3 =30) [211 ■ In these cases tan/3 could also be measured 
to a precision ~ 5% [SEl- At the largest values of mA-, where only the A^ -^ tt is 
available, the accuracy on tua will be typically of the order of a few %, limited by 
energy calibration systematics. 



3.2.4- Other Scenarios Measurements of SUSY parameters and hence dark matter 
properties in other non-CMSSM scenarios have been studied in far less detail. If the 
LSP is a wino then evidence could be provided by studying the branching ratios of the 
X2 and x5) as well as measuring the full neutralino mass matrix. Given a plausible model 
such as that studied in Ref. ^j it may then be possible to measure sufficient parameters 
to estimate Q^h"^. Similarly if the LSP is a sneutrino then this will significantly alter 
the cascade decay chains used in the CMSSM analysis but not prevent observation of 
the kinematic end-points used to constrain masses or model parameters. If a GMSB 
model is realised in nature then the prospects for direct detection (Sectional are slim, 
but studies at the LHC will be fruitful with a wide variety of signatures and mass 
constraints accessible jUj. This should enable the GMSB model parameters to be 
estimated and, if the lifetime of the NLSP can be measured, it should be possible to 
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estimate the SUSY breaking scale and hence the mass of the gravitino LSP. 

4. Cold Dark Matter and Extra Dimensions 

Extra dimensions offer an alternative solution to the hierarchy problem. Instead of 
introducing a cancellation of the divergences as in supersymmetry, extra dimensions 
bring the effective unification scale from the Planck to the TeV scale. Of the several 
models which have been proposed, Universal Extra Dimensions [HI] are of particular 
interest here since they offer a cold dark matter candidate. The simplest realisation of 
UED has all the SM particles propagating in a single extra dimension of size i?, which is 
compactified on a Si/ Z2 orbifold. This single parameter determines the phenomenology 
of the model. Each Standard Model particle has an infinite tower of Kaluza Klein 
(KK) excitations, with each level one excitation possessing a mass equal to 1/ R at tree 
level. This mass degeneracy is broken only by radiative corrections. A further specific 
feature of UED is the conservation of the Kaluza-Klein number at tree level. Radiative 
effects innHnHESj break KK number conservation down to a discrete conserved quantity 
known as KK parity, given by (—1)" where n is the KK level. The level one KK states 
in UED models with conserved KK parity therefore strongly resemble the particles 
appearing in SUSY models with conserved R-parity and indeed it can be difficult to 
distinguish such models experimentally. The lightest KK partner (LKP) at level one is 
often neutral and has negative KK parity causing it to be stable on cosmological time 
scales. The LKP is therefore a plausible dark matter candidate again resembling LSP 
candidates. 

In minimal UED models, the LKP is the KK partner of the electroweak hypercharge 
gauge boson B^ ^ denoted by B^^' jHSj- The relic density of LKP dark matter lies in the 
range favoured by WMAP data if 1/ R ^ 0(1 TeV) [02] • The z/^^^^ state, another possible 
candidate, is already excluded by present limits from direct detection experiments |67j . 
In the early universe B^^^B^^^ annihilation proceeds through the s-channel and is more 
efficient than in the SUSY neutralino case. The cosmologically interesting mass scale 
for the LKP is therefore significantly larger than that for supersymmetric candidates 
(see Figure ITUI). although it still lies well within the anticipated LHC reach of ~ 1.5 TeV 
|68j . Co-annihilation with the e)j state is possible if its mass is almost degenerate with 
that of the LKP however the UED co-annihilation cross-section is not sufficiently large 
to reduce significantly the 5*-^^ density. The net result is that the dark matter density 
increases in the case of nearly degenerate states, since the B^^^ and the e\^ decouple at 
about the same freeze-out temperature. 

Kaluza-Klein dark matter also offers interesting prospects for direct or indirect 
detection at space- and ground-based experiments [HZl IMj- Direct detection (see 
Section El for the SUSY case) proceeds through KK quark and Higgs exchange between 
relic LKPs and atomic nuclei and can give sensitivity to LKP masses in excess of 1 TeV 
for tonne scale experiments and favourably small values of the level one quark and B^^^ 
mass difference. 
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Figure 10. Relic darli matter density as a function of the UED mass scale rriKK = 1/-R 
corresponding to the mass of the level one KK states. The WMAP result is shown by 
the horizontal band. The solid line includes only the B^^' contribution. The dashed 
lines also include the effect of co-annihilations with one and three generations of Icpton 
excitations for mass splittings of 1 % and 5 % (adapted from Ref. [M|)- 



There are three main indirect signatures which can be sought. The first consists 
of a peak of mono-energetic positrons which could be produced in B^^^B^^'^ -^ e~^e~ 
annihilation. The AMS experiment on the ISS should be able to distinguish positrons 
from photons up to about 1 TeV and consequently could provide a means to detect 
LKP dark matter for masses up to ~1000 GeV. A second signature is represented 
by an excess of photons, with a continuum spectrum, which could probe mass scales 
up to about 600 GeV. Finally high energy neutrinos could be detectable at foreseen 
experiments, with a sensitivity up to about 400 GeV. 

Non-supersymmetric grand unified theories with warped extra dimensions also 
introduce a stable Kaluza-Klein state, which could be a candidate for cold dark matter. 
Such a state arises from imposing baryon number symmetry which is spontaneously 
broken at the Planck scale [70j . It has the quantum numbers of a right-handed neutrino 
and cannot decay into SM particles. Depending on the other parameters of the model, 
there is a vast range of masses for this particle to produce the required dark matter 
density, from a few tens of GeV up to about 1 TeV or higher. At the LHC, a possible 
signature of such a scenario is provided by the production of the next-to-lightest KK 
particles which would have a sufiiently long lifetime to travel through the detector 
volume and be recorded as heavy stable particles. 

The main feature of UED models, as well as SUSY models, which leads to the 
existence of a plausible dark matter candidate, is the existence of a conserved quantum 
number distinguishing the new states from those appearing in the SM. Conservation 
of these quantum numbers in general provides these models with additional attractive 
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features jTTj. This idea has been generahsed to solve the so-called "little" hierarchy 
problem suffered by many models predicting new physics at the TeV scale jTT] . In doing 
so a candidate dark matter particle is invariably produced. One test case considered 
by the authors of Ref. [7T] was that of a little Higgs model, in which a Z2 symmetry 
referred to as "T-parity" is imposed. The result is a dark matter candidate which is 
either the new B' gauge boson or a new SU{2)w singlet or triplet scalar particle. The 
annihilation of the former to SM particles is not chirally suppressed, unlike in e.g. SUSY 
models, and so indirect detection signatures are expected to be enhanced. Singlet scalar 
candidates have very small couplings to other states and must therefore be of very low 
mass (< 100 GeV) to give an acceptable relic density i7lj. States with a significant 
triplet component can interact with light gauge bosons, providing a mechanism for their 
annihilation in the early universe, and consequently their cosmologically allowed mass 
range is expected to be higher (> 500 GeV) ji72j . 

5. Direct Detection of DM particles 

5.1. Motivation 

Although LHC experiments can provide supporting evidence for the existence of cold 
dark matter particles, their limited sensitivity to the lifetime of neutral particles (e.g. 
r < 1 ms for Xi ~^ G'~f in GMSB SUSY models in ATLAS [H]) prevents them from 
providing conclusive proof. Similarly it is very difficult for astroparticle experiments to 
prove that any observed signal is due to a particular candidate, such as the LSP, rather 
than some other exotic form of new physics {e.g. UED). This means that in order to 
obtain a discovery of a specific new dark matter candidate consistent signals will need 
to be observed at both colliders such as the LHC and in astroparticle experiments. 

There are two classes of astroparticle searches for cold dark matter; indirect searches 
seeking evidence of dark matter annihilation products such as neutrinos generated in 
the earth, sun or galactic center, or gamma rays or anti-matter generated in the galactic 
halo (see e.g. Ref. ^| and references contained therein), and direct searches sensitive 
to the elastic scattering of cold dark matter particles from atomic nuclei in earth-bound 
detectors [73]. In this review we concentrate on the latter, however it should be noted 
that the former can also provide very strong constraints and indeed may already have 
provided first evidence for particle dark matter in the form of an excess of positrons 
from the galactic halo 



5.2. Basics of Direct Detection 

In direct searches, evidence for the existence of cold dark matter particles is provided 
by an observation of an anomalous source of recoiling nuclei inside a detector. The 
precise form of the recoil energy spectrum of the nuclei depends on factors arising 
from astrophysics (the local density and velocity distribution of incident dark matter 
particles) , nuclear physics (form- factors and coupling enhancement factors for scattering 



Cold Dark Matter and the LHC 25 

at finite momentum transfer q) and particle pliysics (tlie particle mass and interaction 
cross-section). In general, the spectrum is expected to fall quickly with increasing energy 
due to the small average kinetic energy of the population of dark matter particles trapped 
within our galaxy (the galactic "halo" ) . This leads to a need for detectors with extremely 
low (~ keV) nuclear recoil energy thresholds. In addition, the nuclear scattering cross- 
section is strongly enhanced by a factor oc A^ (nuclear target mass squared) for spin- 
independent interactions (see below), favoring target materials incorporating high mass 
(large Am) nuclei. 

The main background to signal events due to nuclear recoils in these detectors 
arises from electron recoils generated by Compton scattering of gamma radiation from 
naturally occurring radio-nuclides in, and around, the target. It is therefore essential 
that detectors be shielded from external sources of gamma radiation, and that they 
be constructed from materials of the highest possible radio-purity. It is also most 
advantageous if detectors are designed to discriminate between nuclear recoil signal 
and electron recoil background events. 

In addition to the electron recoil background there is also a potential source of 
irreducible background events from elastic scattering of neutrons from nuclear spallation 
by cosmic ray muons and decay of radio-nuclides from natural U and Th chains in the 
detector and its environment. The flux of spallation neutrons can be reduced to very 
low levels by installing detectors deep underground where the muon flux is minimal. 
Neutrons from radioactive decay are generally of lower energy and can be absorbed 
using hydrogenous shielding materials or tagged using active neutron vetoes. 

5.3. Direct Searches and SUSY 

Direct searches for cold dark matter rely upon the existence of weak-scale couplings 
between dark matter particles and quarks and gluons in the nucleus. The nature of 
these couplings depends strongly on the identity of the dark matter particle with scalar, 
axial- vector and vector couplings being important for different candidates |[2^. Axial- 
vector currents couple to the distribution of spin in the nucleus and hence contribute to 
the spin- dependent elastic scattering cross-section. Scalar currents by contrast couple 
effectively to the mass distribution of the nucleus and contribute to the spin-independent 
elastic scattering cross-section. For non-Majorana dark matter particles vector couplings 
tracing the charge distribution of valence quarks in the nucleus can further contribute 
to the spin-independent cross-section j73]. It is important to note that the spin- 
independent cross-section is a measured quantity and is, in this case, not equivalent 
to the scalar cross-section, which must be calculated using model assumptions. 

In gravity-mediated MSSM cold dark matter models the LSP can be the lightest 
neutralino, a sneutrino, or possibly a light gluino (see Section EH) . As noted above, the 
last possibility is heavily constrained and will not be discussed further. The sneutrino 
has a large coupling to ordinary matter and hence is almost excluded jSZ] by direct 
searches due to a large predicted vector coupling contribution to the spin-independent 
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scattering cross-section. This leaves the Majorana Xi ^^ the most favored MSSM cold 
dark matter candidate. 

In CMSSM models when the LSP is neutral it is always the Xi- At tree- level 
the scalar x5"iiucleus interaction receives contributions from q — Xi couplings involving 
t-channel Higgs exchange or s-channel q exchange. At 1-loop further significant 
contributions can arise from diagrams coupling the Xi to a gluon [76^. Axial- vector 
interactions arise from t-channel Z° exchange and s-channel q exchange. The magnitudes 
of the scalar and axial-vector cross-sections depend on a number of SUSY parameters 
including the Xi mass and composition {i.e. Mi, tan/5 and /i), the squark masses and 
the heavy Higgs mass niA- In order to obtain a model-independent estimate at the LHC 
of the expected event rates in direct search experiments it will therefore be necessary to 
perform a wide range of SUSY particle mass measurements. 

Within the CMMSM, typical values for the scalar x^'i^ucleon interaction cross- 
section range from ~ 10~^ pb to ~ 10~^ pb in the focus point region, where XiX?^ 
couplings are enhanced by the Higgsino content of the Xi and Aqq couplings by the 
large value of tan/3 [77j, to < 10^^° pb in the co-annihilation tail [2H1- When yu < it 
is possible for the u and d quark contributions to the scattering to interfere destructively 
reducing the cross-section still further fIE\ I7H] . Large values of the scalar cross-section 
are also obtained in regions with small rriQ and mi/2 where the squark mass is small 
leading to efficient s-channel scattering. Similar arguments apply to the axial-vector 
x'j'-nucleon cross-section, with values ranging from 10~^ pb in the focus point region, 
where the Higgsino content of the Xi enhances the Z^XiXi coupling, down to 10~^ pb 
in the stau co-annihilation tail, where rrigj^ is large [75] . 

5.4- Prospects for Direct Searches 

The current generation of direct search experiments have attained a sensitivity to the 
normalized spin-independent Xi-nucleon cross section (Tsj ~ 4x 10"^ pb [7^] (Figure ITTj). 
The DAMA experiment operating at the Gran Sasso underground laboratory ^^ 
has claimed evidence for an annual modulation signature consistent with a cold dark 
matter interpretation (m^^o ~ 52 GeV, agi ~ 7.2 x 10^^ pb |H11), however the region 
of parameter space allowed by this data has now been largely excluded by other 
experiments [73 IMl EH IH3| ■ The region allowed by the DAMA data can be extended 
by varying the parameters of the velocity distribution of cold dark matter particles 
trapped within our galaxy [M], or by assuming contributions to the signal from both 
spin-dependent and spin-independent couphngs |HS]- Nevertheless it is not clear that 
these new regions of parameter space are not also excluded by other experiments when 
similar assumptions are made. Current limits on the spin-dependent Xi-nucleon cross- 
section asd assuming a pure Higgsino LSP candidate (see below) are ~ 1 pb [HE]- 

Within the next two years a new generation of high mass experiments capable 
of enhanced electron-recoil background discrimination will come on-line. These are 
sensitive to spin-independent interaction cross- sect ions ~ 10"^ pb after a few years 
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Figure 11. Sensitivities of some running and planned direct detection dark matter 
experiments to the spin-independent elastic scattering cross-section. Full curves 
correspond to limits from existing experiments (ZEPLIN-1 [HS], EDELWEISS f^ . 
CDMS-II (2004) [ZH]), dashed curves to predicted sensitivities of future experiments 
(next generation: CRESST-II W, ZEPLIN-2 [M). CDMS-II EOl and EDELWEISS-2 
[SI; tonne-scale: GENIUS ,97,, XENON (HSI and ZEPLIN-4/MAX [SliS,). The full 
dark (red) region corresponds to the 3 a allowed region from the DAMA experiment 
|84|. The full light regions correspond to predictions by Chattopadhyay et al. |5S) 
and Baer et al. |2H1 in the light of WMAP data. The crosses correspond to neutralino 
masses and cross-sections predicted for post-LEP benchmark CMSSM models |5(J|. 
Adapted from Ref. '8T. 



of running (HHl EOl 1^ 1121 ■ In the longer term, specific proposals are being made 
to build still larger tonne-scale experiments jHSl HH HH UHl HZI capable of probing 
asi ~ 10"^° pb at the lower limit of the CMSSM favored cross-section region (at least 
for /i > 0). These experiments are likely to be particularly challenging because of the 
difficulty of ensuring that neutron-induced background events do not systematically 
limit detector sensitivity. It is likely that the neutron flux from external sources can 
be reduced to acceptable levels with active and passive shielding [OSj. Reduction of 
the neutron flux from U/Th contamination of detector components will require careful 
design and material assay. This background will probably deflne the ultimate sensitivity 
of direct search experiments since even the purest detector construction materials will 
always contain some U/Th chain nuclides [HH]- Improvements in the sensitivity to 
spin-dependent interactions can be expected to be similar to those for spin-independent 



Cold Dark Matter and the LHC 28 

interactions, however the lack of a cross-section enhancement due to nuclear coherence 
[75] means that tonne-scale detectors will only be sensitive to agd ~ 10~^ pb, which is 
only just inside the CMSSM predicted range. 

By 2011, when the LHC should have acquired 100 fb~^ of data, tonne-scale direct 
search experiments should have been running for a few years. By this stage therefore cold 
dark matter models predicting agi ^ 10~^ - 10~^° pb should have either been excluded 
or produced a signal in these experiments. This opens up the exciting prospect of 
complementary studies of particle dark matter using direct search data and exclusive 
measurements at the LHC. 

6. Complementarity of Measurements 

Direct search experiments are complementary to the LHC in several important ways. 
In terms of discovery of SUSY it is clear that direct searches are particularly sensitive 
to focus point models in which the x? acquires a significant Higgsino component. These 
are the same models however which the LHC will find hardest to discover due to the 
high masses and hence reduced production cross-sections of both squarks and gluinos. 
As a result the LHC and direct search discovery contours are roughly orthogonal in 
this region of parameter space and the latter could conceivably find evidence for BSM 
physics at mass scales in excess of 3 TeV with agi as high as 10~^ pb J2H]- By contrast, 
in the co- annihilation tail the sensitivity of direct search experiments drops off rapidly 
with rising mi/2 as rrigj^ increases, and here the LHC stands a much better chance of 
making a discovery. 

In the event of discovery of signals at both the LHC and direct detection 
experiments further valuable cross-fertilization is possible. If the LHC experiments 
can measure the SUSY model parameters {e.g. mo, rni/2 etc. in the CMSSM) 
using the techniques described in Section 13.21 then a model-dependent estimate of 
the spin-independent and spin- dependent elastic scattering cross-sections can also be 
obtained. It may also be possible to measure sufficient weak-scale SUSY parameters, 
including the neutralino mass matrix parameters together with rriA and rrigj^, to obtain 
a "model-independent" estimate of the cross-sections for comparison with direct search 
measurements. Note that the precision of these measurements will itself be ultimately 
limited by astrophysical model-dependence due to the unknown distribution of dark 
matter within the galactic halo. It is important to remember however that the direct 
search parameter space is defined both by the interaction cross-section and the dark 
matter particle mass. The kinematics of the elastic scattering process mean that for x^ 
masses < Am (the mass of the target nucleus) the shape of the energy spectrum of signal 
nuclear recoils for scalar interactions is strongly correlated with the mass of the Xi (fo'^ 
larger x1 masses the energy spectrum tends to a limiting distribution insensitive to the 
mass value). This correlation is furthermore rather model- independent. Consequently 
for x^ masses < 180 GeV (~ the mass of a tungsten nucleus, the heaviest nucleus used 
in current or proposed experiments [02]) it should be possible to measure the Xi mass 
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rather accurately without recourse to any particular model framework. This will then 
allow direct comparison with LHC results. 

It is interesting to note that a good case study of a direct search mass measurement 
already exists in the form of the analysis performed by DAMA of their claimed annual 
modulation signal jHUl IHH IHSI • The soft energy spectrum observed in this experiment 
is claimed to be consistent with a small Xi mass considerably less than that of the ^^^I 
target nuclei. Under this assumption therefore an accurate measurement of the Xi mass 
can be performed giving a precision ~ 20 % {52^1^ GeV ^j). 

LHC-direct search complementarity is likely to be even more pronounced if a spin- 
dependent dark matter signal is observed. In this case, there is a strong dependence of 
the relative magnitude of the axial- vector couplings of the Xi composition to protons and 
neutrons [011 ^^^^ ^^^ ^^^^ ^^ significant model-dependence in the overall Xi-nucleus 
(or normalized Xi-nucleon) cross-section. Only in the pure Higgsino case is this model- 
dependence lifted. Direct search experiments can absorb this model-dependency into 
a redefinition of measured cross-sections in terms of cross-sections for scattering from 
protons and neutrons separately 1^3 ■ Only with the aid of LHC measurements of e.g. 
the neutralino mass matrix however can these results be combined to obtain an overall 
experimental measurement of the spin-dependent scattering cross-section. 

The variation of the relative magnitude of Xi-P and Xi'T^ axial- vector couplings with 
Xi composition also leads to some residual model dependency in the shape of the energy 
spectrum obtained from spin-dependent scattering. The form-factor S{q) relating x^- 
nucleus scattering at finite momentum transfer q to that at zero momentum transfer 
contains terms depending on isoscalar and isovector parameters Oq and Oi related to 
Xi-p and Xi-n coupling parameters ap and a„ |innj : 

S{q) = alSooiq) + ajSuiq) + aoai5'oi(g) (2) 

where 

ao = ap + an, ai = ap — an. (3) 

Here the Sij{q) are independent constituent nuclear form-factors. The differing q 
dependence of these Sij{q) form-factors means that S{q) differs for differing ap/an and 
hence differing x^ composition. This leads to a SUSY model dependency in the shape of 
the overall energy spectrum. This model-dependency is generally rather small (< 10 % 
for low energy recoils in e.g. ^^^I) however the accuracy of direct search measurements 
of asd and m^o for comparison with LHC predictions could conceivably be improved 
with an estimate from LHC data of the neutralino mass matrix. 

7. Conclusions 

Recent astrophysical measurements from WMAP and other experiments have 
strengthened the case for cold dark matter in the universe. Over the next decade the 
hypothesis of particle cold dark matter will be stringently tested in a number of different 
ways both with colliders such as the LHC and with non-accelerator experiments such as 
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direct dark matter searches. If SUSY or other BSM physics is discovered at the LHC then 
measurements of the properties of the new particles can be used to predict the properties 
of dark matter signals in CMB and astroparticle experiments. Similarly if dark matter is 
discovered in a direct search experiment model-independent measurements of quantities 
such as the mass of the dark matter particle can be compared with signals observed 
at the LHC. Only by combining these results can we confirm both the existence of 
non-baryonic cold dark matter and its identity. 
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